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From  the  start,  my  objective  was  to  develop  a low-cost 
method  of  increasing  the  sampling  rate  of  a rate  aided 
pointing  and  tracking  system  by  using  information  already 
available  within  the  system.  After  determining  the  best 
of  alternative  pseudo-measurement  schemes,  the  system  proved 
to  be  very  effective  against  highly  maneuverable  targets. 
Although  only  a two  dimensional  analysis  with  simplified 
dynamics  was  used,  it  is  expected  that  the  general  results 
will  carry  over  in  a satisfactory  manner  to  three  dimensional 
systems  with  more  complicated  dynamics. 

My  sincere  appreciation  goes  to  Prof.  James  E.  Negro 
for  his  advice  and  patient  guidance  in  this  effort.  Most 
of  all,  I must  express  my  gratitude  to  my  wife,  Berni,  for 
her  encouragement  and  endurance  of  many  lonely  hours  during 
the  preparation  of  this  work. 
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Abstract 

This  thesis  investigates  the  use  of  Internally  generated 
pseudo-range  measurements  to  increase  the  sampling  speed  of 
a rate  aided  pointing  and  tracking  system  known  as  Aided 
Track.  Six  separate  methods  for  generating  pseudo-range 
measurements  were  developed  for  a two  dimensional  analysis 
of  the  problem.  Five  target  scenarios  of  various  complexity 
were  used  to  test  the  systems.  The  results  obtained  indicate 
that  the  use  of  a pseudo-range  measurement,  generated  by  using 
the  internal  estimate  of  range  provided  the  best  performance 
of  all  methods  tested.  This  method  demonstrated  a definite 
improvement  over  the  slower  system  when  tracking  highly 
maneuverable  targets. 


PSEUDO- RANGE  MEASUREMENTS 
USED  TO  INCREASE  THE  SAMPLING 
RATE  OF  THE  AIDED  TRACK  ALGORITHM 


I . Introduction 

Background 

The  problem  of  accurately  pointing  a gimbaled  sensor 
at  a target  and  then  precisely  tracking  that  target  has 
always  been  important.  Conventional  methods  of  pointing 
and  tracking  have  an  inherent  pointing  error  caused  by 
dynamic  lag.  This  error  is  unacceptable  for  many  applica 
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Tracker 
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Fig.  1.  Pointing  and  Tracking  System  Model 

A new  approach  to  the  pointing  and  tracking  problem 
was  developed  by  Fitts  (Ref  1,2,3).  In  this  approach, 


known  as  Aided  Track,  an  auxiliary  signal  Waid  is  input  to 
the  rate  command  signal  as  shown  in  Figure  1,  Ideally,  if 
is  chosen  properly,  the  tracking  error  e would  be 
forced  to  zero.  The  Aided  Track  concept  is  covered  in 
greater  depth  in  Chapter  II. 

In  order  to  generate  Aided  Track  uses  recursive 

Kalman  filtering  techniques  to  weight  and  incorporate  the 
measurements  of  range  and  gimbal  rotation  rates.  Nominally, 
the  gimbal  rate  information  can  be  provided  at  0.01  second 
intervals  by  gyros  attached  directly  to  the  girnbaled  sen- 
sor. Whereas,  accurate  range  information  is  only  available 
every  0.1  seconds  from  a laser  rangefinder.  Aided  Track, 
as  developed  by  Fitts,  operates  at  the  slower  10  samples/ 
second  data  rate  due  to  the  laser  rangefinder  limitation 
(Ref  1:37). 


Purpose  of  the  Study 

At  the  current  sample  rate,  nine  out  of  ten  measure- 
ments of  gimbal  rate  are  being  ignored.  Logically,  it 
would  seem  that  if  these  nine  measurements  could  be  incor- 
porated in  some  way,  that  an  improvement  in  the  system 
could  be  realized.  The  simplest  way  to  do  this  would  be 
to  increase  the  sample  rate  of  the  laser  rangefinder,  how- 
ever, this  is  not  possible  due  to  heat  dissipation  limita- 
tions. Therefore,  some  other  means  of  obtaining  range 
measurements  to  couple  with  the  nine  gimbal  rate  measurements 
must  be  used.  The  purpose  of  this  study  is  to  determine  the 
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best  method  of  internally  generating  pseudo-range  measure- 
ments to  combine  with  the  extra  gimbal  rate  measurements 
and  to  determine  if  and  how  much  of  an  improvement  in  track- 
ing performance  can  be  gained. 

Out  line  of  the  Report 

In  this  study,  a two  dimensional  analysis  of  the  prob- 
lem is  presented.  This  approach  allows  for  a more  simpli- 
fied computer  simulation  with  the  results  being  easily 
related  to  the  more  complicated  three  dimensional  case. 

The  actual  computer  simulation  developed  is  listed  in 
Appendix  A. 

Chapter  II  presents  a discussion  of  Aided  Track  and 
some  of  the  rationale  behind  it.  The  actual  transforma- 
tions and  associated  equations  are  presented  in  Chapter 

III. 

Chapter  IV  then  describes  each  modification  of  the 
program  that  was  used  to  generate  the  pseudo-range  measure- 
ments. The  rationale  for  each  modification  is  also 
presented  along  with  the  problems  encountered  and  how  they 
were  remedied. 

Next,  Chapter  V presents  an  analysis  of  the  results 
which  shows  an  agreement,  in  part,  with  those  obtained  by 
Fitts  (Ref  1,2,3).  Any  attainable  improvement  in  tracking 
performance  is  shown  to  be  dependent  on  target  motion. 

Finally,  Chapter  VI  presents  the  conclusions  and 
recommendations  regarding  the  overall  analysis. 
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II.  Aided  Track 


The  analysis  of  Aided  Track  presented  here  is  for  a 
simplified  two  dimensional  model.  Figure  2 depicts  a 
model  of  the  azimuth  channel  of  a pointing  and  tracking 
system  wherein  the  stabilization  loop  is  modeled  by  the 
transfer  function 
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which  is  a second  order  approximation  of  the  actual  stabi 
1 i za  ti on  1 oop . 
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Fig.  2.  Azimuth  channel  of  Pointing  and 
Tracking  System  (Ref  4:41) 

I 

The  natural  frequency  wn  and  the  damping  ratio  c of  the 

n 

stabilization  loop  are  220  radians/second  and  0.7  re- 
spectively (Ref  2:40). 
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The  tracker  and  proportional  plus  inteqral  compensa- 
tion is  accurately  modeled  by  the  transfer  function 

G t ( s ) = ¥S/m  * (2) 


For  this  analysis,  the  tracker  sensor  sampling  period  (T) 

_ ? 

is  0.01  seconds.  Therefore,  a gain  Ky  of  700  seconds 
and  a frequency  Wy  of  18  radians/second  give  suitable 
closed  loop  characteristics  (Ref  3:41-44). 

Ideally,  one  would  like  W ^ to  be  chosen  so  that  the 
tracking  error  e is  forced  to  zero.  From  inspection  of 
Figure  2,  one  sees  that  zero  error  implies 

Waid  (S)  = SnT (S)/GS(S)  (3) 

Substituting  for  Gs(s)  yields 

Waid  (S)  3 [S  + (2c/Wn)S2  + S3/Wn2]nT(S)  (4) 

Now,  taking  the  inverse  transform  of  both  sides  of  Eq 
(4)  yields 

# •• 

Waid(t)  = ny  + 2r,  n y / p + higher  order  terms  (5) 

The  problem  of  generating  ^ ^ is  that  the  derivatives 

of  the  target  angle  specified  in  Eq  (5)  are  not  available 
from  direct  measurements.  However,  Fitts  determined  that 
sufficiently  accurate  estimates  of  these  derivatives  could 
be  generated  by  using  a Kalman  filter  external  to  the  servo 
loop  (Ref  1).  Eq  (5)  then  becomes 
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(6) 


Waid  (t)  = ~T  + C,A 

where 

Ga  = 2r,/Wn  + 3.5  ATT/2*  (7) 

is  the  quadratic  stabilization  loop  correction  term.  The 
first  term  of  the  qain  follows  from  Eq  (5).  The  second 
term,  however,  was  selected  by  Fitts  to  compenstate  for 
the  lag  due  to  the  discretization  of  the  actual  aided 
track  signal  (Ref  1:16). 

In  order  to  generate  the  estimates  of  gimbal  rate  and 
acceleration,  Fitts  elected  to  use  a linear  Kalman  filter. 
This  approach  requires  several  transformations.  First,  the 
measurements  Rm,  nm,  nm  along  with  the  estimate  of  range 
rate  from  the  filter  must  be  transformed  into  an  inertial 
coordinate  frame.  The  measurement  sample  interval  for  the 
Fitts  model  is  AT  = 0.1  seconds.  This  interval  will  be 
decreased  in  this  study  to  AT  = 0.01  seconds  to  incorporate 
pseudo-range  measurements.  The  filter,  however,  operates 
at  ATT  = 0.01  seconds  already,  in  order  to  provide  a 
smooth  Waid  command  signal  compatible  with  the  servo  loop 
8.0  Hz  frequency.  Estimates  of  the  target  position, 
velocity,  and  acceleration  in  an  inertial  frame  are  output 
from  the  filter.  These  estimates  must  then  be  transformed 
back  into  polar  coordinates  to  provide  the  estimate  of 
range  rate  for  the  filter  and  estimates  of  gimbal  rate  and 
acceleration  for  generation  of  1^.^. 
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This  chapter  has  presented  a very  basic  description  of 
the  Aided  Track  concept  as  developed  by  Fitts.  It  again 
should  be  emphasized  that  the  design  rationale  is  the 
desire  to  use  a linear  Kalman  filter  with  precomputed 
gains  rather  than  on-line  gain  calculations  required  in 
an  Extended  Kalman  filter  approach  for  non-linear  systems. 
Details  of  the  Aided  Track  approach  are  included  in  Chapter 
III. 


III.  Si mulation  Model 

The  Aided  Track  model,  as  discussed  in  Chapter  II,  is 
depicted  in  Figure  3.  In  this  chapter,  the  reference 
frames  in  which  the  target  motion  is  modeled,  are  defined. 
Also,  the  equations  representing  the  truth  model  and  Kalman 
filter  are  presented  along  with  their  associated  transfor- 
mations. 

Ref e re nee  F rames 

Fitts  defines  several  reference  frames  in  order  to 
accommodate  tracker  motion  in  his  three  dimensional  analy- 
sis (Ref  1,2,3).  For  purposes  of  this  study,  the  tracker 
will  be  assumed  non-translating  with  respect  to  inertial 
space,  thereby  requiring  only  two  coordinate  frames  to  be 
defined.  The  orientation  of  these  reference  frames  is 
arbitrary,  but  usually  chosen  for  convenience. 

An  inertial  coordinate  system  Yj,  is  defined  to  be 
fixed  on  the  surface  of  the  earth  with  the  tracker  located 
at  the  origin.  For  convenience,  the  Y-j  axis  is  defined  to 
be  pointing  directly  at  the  target  at  t = 0. 

A similar  inertially  non-rotating  coordinate  system 
X1 , is  defined  for  the  Kalman  filter.  However,  it  is 
initialized  when  the  filter  is  turned  on  so  that  the 
axis  is  pointing  directly  at  the  estimated  position  of 
the  target.  The  target  geometry  can  now  be  described  as 
illustrated  in  Figure  4.  The  angle  n0  is  fixed  at  the 
Instant  the  filter  is  activated  and  is  used  to  make  trans- 
formations from  the  filter  frame  back  to  the  truth  model 
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Figure  3:  Aided  Track  Model 


inertial  frame.  This  angle  is  used  only  as  a simulation 
artifice  to  relate  the  filter  frame  to  the  simulation  frame. 
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Fig.  4.  Target  Geometry 


Truth  Model 

The  truth  model  consists  of  the  differential  equa- 
tions that  model  both  the  target  motion  and  the  gimbal 
dynamics.  The  matrix  expression  for  the  truth  model  is 
shown  in  Eq  (3).  The  formulation  of  this  expression  is 


10 


Equation 


developed  in  Appendix  B.  This  matrix  expression  of  the 
truth  model  is  solved  by  a differential  equation  solving 
routine  using  a variable  step  size  down  to  1.0  E-6  seconds. 
Output  is  actually  provided  every  hundreth  of  a second  in 
order  to  simulate  the  100  Hz  sample  rate  of  the  servo 
loop.  The  laser  rangefinder  is  then  simulated  by  sampling 
the  true  target  position  every  0.1  seconds.  Similarly, 
the  gimbal  rate  measurement  is  simulated  by  sampling  W ^ 
every  0.01  seconds.  Transformation  T«  in  Figure  3 also 
uses  W j in  the  equation 


where 


nm(K+l) 

K+l 


* *TT  Wcmd 
ATT  + tk 


(9) 


to  generate  nm-  At  filter  turn  on  time  tQ>  the  angle  n()1 
is  initialized  to  zero. 


Transformation  T-| 

In  order  to  utilize  a linear  Kalman  filter  as  men- 
tioned in  Chapter  II,  the  polar  coordinate  measurements 

• • 

(Rm,  R , nm,  n|n)  must  be  transformed  into  a planar  coordi- 

• • 

nate  pseudo-measurements  (X^m,  X-j n) , X2|r|,  X2m).  The 
measurement  of  range  rate  Rn),  and  target  angle  nm»  are 
not  readily  available.  However,  the  target  angle  n|n  is 
provided  as  was  shown  in  Eq  (9),  and  range  rate  Rm 
initally  can  be  obtained  from  a weighted  combination  of 
the  last  three  laser  range  measurements  by  the  equation 


whore  t is  the  filter  turn  on  time  and  AT.  = 0.1  seconds 
o L 

(Ref  1:39).  Subsequently,  an  estimate  of  range  rate  ft 
is  available  from  the  filter.  The  pseudo-measurement 
generation  can  now  be  expressed  as  follows: 


1 m 

= Rm 

COS 

M.n 

2 in 

* Rm 

SIN 

nm 

1 m 

= ft 

m 

COS 

nm 

Rm 

WCMD 

SIN  n 

m 

2m 

= ft 

m 

SIN 

nm  + 

Rm 

WCM0 

COS  n 

m 

(ID 


The  pseudo-measurements  are  now  ready  to  be  processed  by 
the  filter. 


Kajima  n F i Iter  f guations 

The  Kalman  filter  propagates  target  position,  velocity, 
and  acceleration  estimates  forward  in  time  up  to  the  next 
measurement  update  time.  At  that  point,  the  estimates  are 
optimally  combined  with  the  measurements  to  provide  the 
new  best  estimate  of  target  position,  velocity,  and 
acceleration.  The  propagation  and  update  equations  are 
separate  processes  and  are  presented  accordingly. 

Propagation  fguations.  The  propagation  process 
utilizes  the  state  transition  matrix  (Ref  4:356) 


as  ~ 


The  estimates  arc  propagated  by  the  matrix  equation 


X. (tk+AT) 

'x,ukr 

X(tk+AT)  = 

^(tk+AT) 

-A 

_VV&TL 

_v  v_ 

i = l ,2 

where  the  minus  sign  designates  the  estimates  just  prior 
to  incorporation  of  the  next  measurement  and  the  plus  sign 
designates  the  estimates  just  after  the  last  update. 

Update  Equations.  Before  the  estimates  can  be 
updated,  the  appropriate  Kalman  gains  must  be  computed  by 
the  following  equations  (Ref  2:35): 


P(K+1)'  = A P(K)+AT  + a (3x3  matrix)  (14A) 
G ( K+l  ) = P(K+1 )'HT[HP(K+1 )'HT+R]'1 (3x2  matrix)  (14B) 
P(K+1)+  = [I  - G ( K+ 1 ) H]P(K+1)'  (3x3  matrix)  (14C) 


K+l  = tk  + AT, 

K « 

= 0, 

1,2,  . 

• 

• 

H can  be  referred 

to 

as 

the  measurement  matrix. 

Sine 

position 

and  velocity 

are 

the 

only 

pseudo-measurements,  H 

becomes 

the  2x3  matrix 

“l 

0 

o" 

H ■* 

(15) 
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The  Kalman  gain  matrix  G is  a 3x2  matrix  expressed  as 
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^ is  a 3x3  matrix  describing  the  covariance  of  assumed 
unknown  target  motion  or  assumed  modeling  error.  For  this 
case  where  target  acceleration  is  assumed  constant  over  the 
sampling  interval,  target  motion  can  be  modeled  as  follows: 


X(t)  = F X(t)  + b w(t) 


x(t) 

0 

1 

o“ 

x ( t j 

'o' 

v(t) 

= . 

0 

0 

1 

v (t ) 

+ 

0 

MO 

0 

L 

0 

0 

A(t) 

J. 

w(t) 


(17) 


where  w(t)  is  a zero  mean,  Gaussian,  random  white  noise 
sequence  modeling  the  unknown  change  i n acceleration.  For 
a discrete  time  simulation 


E[w(t)w(s)]  = q«  (t-s)  (1R) 

where  q is  considered  the  strength  of  the  white  noise 
sequence  or  a measure  of  the  uncertainty  in  the  target 
model.  If  the  target  acceleration  is  constant  then  q = 0, 
otherwise  q would  increase  appropriately.  The  convariance 
equation  for  the  target  model  would  be 

P(tk+AT)  = AP(tk)AT  +jATA  b q (t ) bTATdi  (19) 


Approximating  the  integration  process  with  an  Euler  Integra 
tion  for  a discrete  time  system  yields  Eq  (14A)  where 


a - 

A b q 

bT  at  at 

A i /4 

aT3/2 

aT2/2 

s 

aT3/2 

AT2 

AT 

aT2/2 

AT 

1 

q AT 


(20) 


(£  then  is  a function  of  the  sampling  interval  AT  and  q the 
uncertainty  of  the  target  acceleration  profile.  Fitts 
(Ref  4:356-357  ) elected  to  normalize  the  Q,  R_,  and  £ 
matrices  and  ignore  products  of  AT  since  AT  <<1  . There- 
fore C[  reduces  to 

a = 


0 0 0 

0 0 0 

0 0 Q 
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(21) 


This  approximation  is  logical  since  target  acceleration  is 
the  only  state  that  was  assumed  constant  over  the  interval 
Through  experimentation,  Fitts  determined  Qgg  = 1.0  to  be 
optimal  against  maneuvering  targets  when  AT  = 0.1  seconds 
(Ref  2:71  ). 

R is  the  measurement  convariance  matrix.  The  values 
expressed  in  it  are  approximated  by  processing  the  con- 
variance  of  the  laser  rangefinder  and  the  gimbal  rate 
sensor  through  Eq  (11),  to  get 
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1 0 

0 R 


2? 


(22) 


Since  the  velocity  measurements  are  dependent  on  estimated 
range  rate,  the  R ^ term  cannot  be  explicity  calculated. 
Through  experimentation,  R^  = 1C  was  found  to  be  optimal 
(Ref  5 : B 3 5 ) . 

P is  the  covariance  matrix  of  the  estimate  errors. 
Taking  into  account  that  Rn)(tQ)  and  Rm(t0)  are  correlated 
by  Eq  (10),  and  that  the  pseudo-measurements  are  correlated 
indirectly  by  Eq  (11),  Fitts  calculated 


i 

1 .5/ATL 

0" 

£<£0>  - 

1.5/aTl 

f>.5/AT(  2 

0 

(23) 

0 

0 

1_ 

where  a T^  is  the  sampling  period  of  the  laser  rangefinder 
(Ref  4:361 ). 

Eq  (23)  is  derived  using  the  orientation  of  the 
axis  at  filter  turn  on  and  the  definition  of  the  covariance 
of  a zero  mean  process. 


Define 

■Rm  - E<R„J 

v 

■ 

I = 

e 

(24) 

% - E<V1 
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then 


(25) 


Now,  from  Eq  (10),  assuming  Rn1(tQ),  Rni^to"AT^’  and  Rm^to"2AT^ 
are  independent  and  E ( Z 1 2}  = E[(Rm  - E { RmJ  )2]  = 1 (be- 
cause of  the  normalization  with  the  R matrix,  = 1). 

E [ Z1  l2)  = ElZ^J  1.5/ATl  = 1.5/ATl  (26) 

E [ Z 2J  = El  (1.5/AT.)2  Z 2 + (2/ AT.  )2Z- 2 + 

2 LI  LI  (27) 

(1/2ATl)2  Z^J  =(6.5/ATl2)  E t Z12}  = 6.5/ATl2 

Once  the  Kalman  gain  matrix  £ is  calculated,  the 
estimates  and  measurements  are  combined  by  the  formula 
(Ref  1:33-34) 

X;i  ( K)+  = X.(K)-  + G(K)[X.m(K)  -Hji(K)"j,  1-1,2  (28) 


Transformation  T2 

A 

Once  the  estimates  are  generated  by  the  filter, 
whether  they  are  updated  or  propagated  estimates,  they  must 
be  transformed  back  into  polar  coordinates  to  generate  the 
Waid  ^or  t^ie  9imba^s  and  ^ for  filter.  The  following 
equations  are  used  in  T2  to  perform  this  transformation. 


I 
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a a p a ? . l /? 

R = (X^  + Xg  )''* 

n * ARCTAN  (X2/X]  ) 

• • • 

A A />  A A a 

R = (X1X1  + X2X2)/R 


A 

n 


A 

n 


A A A /N  A? 

(X 1 X2  - X2X1  ) / R ^ 

tt  ••  m • 

J A A A A A 9 .A  A 

(X1X2  - X2X1)/R"  - 2nR/R 


(29) 


Transformation  T3 

Filter  estimates  of  gimbal  rate  and  acceleration  are 
used  in  to  generate  as  discussed  in  Chapter  II. 

Eqs  (6)  an  (7)  are  restated  here: 


where 


Waid(K)  = A(K)  + Ga  %(K) 

Ga  = 2c/Wn  + (3, 5/?n  ).\TT 


(30) 

(31) 


"aid 

From  a close  look  at  Figure  3,  it  can  be  seen  that 
Waid  1S  dependent  0,1  W d>  which  itself  is  generated  by 
the  sum  of  ^ ^ and  W^.  It  is  important  to  note  that 
there  are  three  possible  methods  of  handling  this  loop  in 
the  system. 

Delayed . A conventional  method  of  handling  it 
Is  to  generate  a new  Wcmj  with  the  Waid  that  was  calculated 
using  the  last  W 1(j.  This  implies  an  inherent  computational 
delay  in  the  system  and  is  the  method  normally  implemented 
in  hardware.  In  equation  form 
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The  sequence  of  events  for  this  case  would  be  as  follows: 

1.  Filter  updates  at  K and  generates  wa-jcj(K)  + 

2.  Truth. Model  integrates  from  K to  K+l  using  W . . 

(K-l)  ' a1d 

3.  Filter  propagates  from  K to  K+l 

4.  Filter  updates  at  K+l  and  generates  W . . 

(K+l )+  a’d 

Instantaneous.  The  second  method  of  handling  this 

system  loop  is  to  assume  a zero  computational  delay.  This 

implies  that  the  filter  equations  could  generate  the 

instantaneously.  This  is  an  ideal  situation  although  not 

practically  implementable.  In  equation  form 

Kcmd<K>  ’ VK)  + "aid'*'*  <33> 

The  sequence  of  events  for  this  case  is  as  follows: 

1.  Filter  updates  at  K and  generates  ^ ^ ( K ) 

2.  Truth  Model  integrates  from  K to  K+l  using 
“aid<K>+ 

3.  Filter  propagates  from  K to  K+l 

4.  Filter  updates  at  K+l  and  generates  W ^ ^ ( K+l  ) + 

Predicted . The  third  method  of  handling  the 
system  loop  is  to  use  the  filter  to  propagate  the  esti- 
mates forward  which  would  generate  a prediction  for  Waid* 
This  method  is  an  attempt  to  solve  the  computational  delay 

and  approach  the  performance  of  the  "ideal" instan- 
taneous method.  In  equation  form 
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“cmd<K'  “ VK)  + “aid<K> 


(34) 


The  sequence  of  events  would  be  as  follows: 

1.  Filter  updates  at  K and  generates  wa-jcj(K)  + 

2.  Filter  propagates  from  K to  K+l  and  generates 
"aid(K+’>- 

3.  Truth  Model  integrates  from  K to  K+l  using 
"aid<K>‘ 

4.  Filter  updates  at  K+l  and  generates  Waf<j(K+l) 

In  this  analysis,  all  three  methods  are  considered, 

The  results  are  presented  in  Chapter  V. 
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I V . Desinn  Improvements 

The  simulation  discussed  in  Chapter  III  must  now  be 
modified  to  produce  the  pseudo-range  measurements  so  as  to 
increase  the  sampling  rate.  Five  basic  modifications  are 
used  along  with  a modification  simulating  actual  range 
measurements  being  available  every  0.01  seconds.  This 
latter  modification  is  performed  so  as  to  determine  the 
optimum  improvement  possible.  This  chapter  will  first 
discuss  what  changes  are  made  in  the  basic  equations  in 
order  to  accept  the  pseudo-range  measurements.  Then,  the 
modifications  are  presented  with  supporting  rationale. 
Finally,  the  scenarios  used  to  test  the  changes  are  dis- 
cussed. 

Filter  Ch an^jes 

Several  of  the  modifications  to  be  mentioned  were 
incorporated  into  the  simulation  with  no  change  to  the 
filter  equations.  Upon  inspection  of  the  output,  a sinu- 
soidal response  was  found  to  be  corrupting  the  system 
response  as  seen  in  Figure  5 (a  plot  of  tracking  error 
versus  time).  A steady-state  perturbation  analysis 
(Appendix  C)  revealed  that  the  inherent  errors  in  the  pseudo- 
range measurements  were  being  amplified  through  the  propa- 
gation process  until  a new  accurate  laser  range  measurement 
could  be  incorporated.  Therefore,  the  filter  had  to  be 
altered  to  make  it  rely  more  heavily  on  the  angle  and 
angular  rate  measurements  and  much  less  on  the  pseudo-range 
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measurements . The  method  chosen  to  do  this  was  to  alter  the 
H matrix  durinq  the  i ncorporat i on  of  the  pseudo-range 
measurements  from 


to 


1 

0 


0 0 

1 0 


H = 


0 0 


0 1 


0 

0 


(35) 


(36) 


The  H.j  i term  was  set  to  zero,  thereby  ignoring  the  pseudo- 
measurements of  target  position  X^,  X2|)i  until  the  next 
laser  range  measurement  could  be  taken.  This  had  to  be 
done  for  all  modifications  where  pseudo-range  information 
was  used.  The  sinusoidal  corruption  was  eliminated  when 
this  change  was  made  to  the  algorithm. 

Another  change  to  the  filter  equations  is  the  rescal- 
ing of  the  Q matrix.  When  the  sampling  rate  is  increased 
the  Q matrix  must  be  scaled  accordingly,  as  shown  by  Eq 
(20).  For  this  study  is  changed  from  1.0  to  0.1. 

A study  of  versus  RMS  pointing  error  performance  was 
accomplished.  The  results  of  this  analysis  are  presented 
in  the  next  chapter. 


Modification  One 

In  modification  one,  range  is  assumed  to  be  constant 
over  the  sample  interval.  This  modification  if  the  simplest 
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to  implement,  and  demonstrates  the  feasibility  of  the  over- 
all study. 

Modification  Two 

Modification  two  propaqates  the  last  laser  range 
measurement  forward  by  use  of  the  internally  generated 
estimate  of  range  rate.  Recall  that  this  estimate  of  range 
rate  is  recomputed  every  0.01  seconds  by  the  filter.  This 
method  was  scrutinized  for  possible  stability  problems  due 
to  the  likelihood  that  errors  in  the  estimate  of  range  rate 
could  propagate  larger  errors.  This  method  did  result  in 
a rather  large  sinusodial  type  response  until  the  change  to 
the  H matrix  was  made.  This  change  eliminated  the  problem. 

Modification  Three 

The  possible  instability  of  the  last  method  was  the 
motivation  for  modification  three.  In  this  method,  range 
is  propagated  by  the  estimated  range  rate  obtained  at  the 
last  laser  range  measurement  time.  By  holding  this  esti- 
mate constant  over  the  0.1  second  interval,  the  propagation 
of  errors  was  minimized.  However  the  sinusoidal  response 
was  still  present  until  the  change  was  again  made  to  the  _H 
matrix. 

Modification  Four 

This  modification  simulates  actual  laser  range  measure- 
ments every  0.01  seconds.  The  motivation  for  this  method 
was  to  determine  the  best  possible  performance  available 


25 


from  increasing  the  sample  rate  and  to  use  this  method  as 
a comparison  for  the  others.  Also,  the  results  from  this 
modification  demonstrate  that  the  sinusoidal  response 
obtained  by  using  the  other  methods  was  due  to  the  propa- 
gation of  errors  and  not  due  to  increasing  the  computation 
rate  of  the  filter  equations  since  this  modification  showed 
no  sinusoidal  corruption  with  the  H matrix  unchanged. 

Mo di f i ca t i o n Five 

The  motivation  for  this  method  again  comes  from  the 
instability  problem  of  modification  two.  In  this  method, 
range  is  propagated  by  a constant  range  rate  obtained  by 
a linear  combination  of  the  last  three  laser  range  measure- 
ments. This  modification  has  a unique  problem  when  noise 
is  introduced  into  the  system.  When  this  is  done,  the 
differenced  range  rate  becomes  more  noisy  than  the  esti- 
mated range  rate.  Again,  the  sinusoidal  corruption  disap- 
peared when  the  H matrix  was  changed. 

Modi fi cat i on  Six 

This  method  of  generating  pseudo-range  measurements 
is  simply  to  use  the  estimated  range  available  from  the 
filter.  Initially,  this  method  was  ignored  because  it  was 
the  method  tried  by  Fitts  (Ref  1:77)  which  reportedly 
demonstrated  poor  performance.  However,  Fitts  considered 
only  flyby  type  scenarios.  Once  the  results  of  Modifica- 


tion Two  demonstrated  the  capabilities  of  this  concept 
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against  maneuvering  targets,  the  use  of  the  estimated  range 
was  again  considered.  This  method  also  required  the  H 
matrix  to  be  changed. 

All  of  the  modifications  mentioned,  as  well  as  the 
altering  of  the  Jj  matrix,  are  handled  by  the  update  equa- 
tions of  the  filter.  Appendix  D contains  a listing  of  these 
changes  as  they  are  implemented  in  the  computer  program. 
Table  I is  a listing  of  the  modifications. 

Table  I 

Listing  of  Modifi ca ti o ns 

Modification Method  used  to  Generate 

One  Range  held  constant  from  last  measure- 

ment. 

Two  Range  propagated  by  estimated  range 

rate . 

Three  Range  propagated  by  constant  range 

rate  obtained  from  last  estimate. 

Four  Laser  range  measurements  every  0.01 

seconds . 

Five  Range  propagated  by  constant  range 

rate  obtained  from  differenced  range 
measurements . 

Six  Estimated  Range. 

Scenarios 

Initially,  a straight  flyby  scenario  was  used  as  shown 
in  Figure  6.  This  scenario  is  very  similar  to  that  used 
by  Fitts  (Ref  1,2,3)  and  was  chosen  so  that  a comparison 
could  be  made  to  insure  that  the  resul ts  agree . Table  II 
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presents  a listing  of  range,  range  rate,  and  angular  rate 
versus  time  for  scenario  one. 

TABLE  II 
SCENARIO  ONE 


Time 
(sec ) 

Range 

(meters  ) 

Range  Rate 
(meters/sec ) 

Angular  Rate 
(radians /sec ) 

Accel erat i on 

( G ~ s ) 

0.  5 

5590 

-984 

.0320 

0 

1.0 

5099 

-981 

.0385 

0 

1.5 

4610 

-976 

.0471 

0 

2.0 

4123 

-970 

. 0588 

0 

2.5 

3640 

-961 

.0755 

0 

3.0 

3162 

-949 

.1000 

0 

3.5 

2693 

-928 

.1379 

0 

4.0 

2236 

-894 

.2000 

0 

4.5 

1803 

-832 

.3077 

0 

5.0 

1414 

-707 

.5000 

0 

5.5 

1118 

-447 

.8000 

0 

6.0 

1000 

-20 

1 .000 

0 

6.5 

1118 

447 

.8000 

0 

7.0 

1414 

707 

.5000 

0 

7.5 

1803 

832 

,3077 

0 

8.0 

2236 

894 

.2000 

0 
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The  second  scenario  used  is  almost  identical  to  the 
first  except  that  at  crossover  the  target  enters  a constant 
radius*  constant  velocity  turn  around  the  tracker  as  shown 


in  Figure  7.  This  scenario  was  chosen  in  order  to  determine 


how  well  this  system  perforins  against  an  unmodeled  accelera- 
tion. The  acceleration  in  the  turn  is  approximately  100  G's 
and  definitely  accentuates  any  difficulties  in  handling 
unmodeled  accelerations. 

These  first  two  scenarios  were  used  to  select  which 
modification  performed  the  best.  Once  this  was  determined, 
the  best  modification  was  tested  against  the  following,  more 
realistic,  scenarios. 

Scenario  three  is  a reduced  version  of  scenario  two, 
with  respect  to  velocity  and  G forces,  as  shown  in  Figure 
8.  Velocity  is  held  constant  at  300  meters/second  and  at 
crossover  the  target  enters  approximately  a 9.2  G turn. 

Table  III  is  a tabulation  of  range,  range  rate,  angular 
rate,  and  acceleration  versus  time  for  scenario  two  and 
three . 

Figure  9 illustrates  scenario  four  which  is  an  adapta- 
tion of  scenario  one.  Velocity  is  held  constant  at  1000 
meters/second  until  t=  4.5  seconds.  At  that  time,  the 
target  undergoes  a deceleration  modeled  by  the  drag  on  a 90 
kilogram  missile  with  a radius  of  15  centimeters  and  a 
coefficient  of  drag  = .25.  The  time  history  of  scenario 
four  is  presented  in  Table  IV. 

Scenario  five  begins  as  in  scenario  one  until  t a 2 
seconds.  At  that  time,  the  target  enters  approximately  a 
16  G turn  toward  the  tracker  until  it  is  pointed  approxi- 
mately at  the  tracker.  The  target  then  rolls  out  and 
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TABLE  III 

SCENARIO  TWO 

Time 

(seconds ) 

Range 
(meters ) 

Range  Rate 
(meters/sec ) 

Angular  Rate 
(radians/sec) 

Acceleration 
( G 1 s J 

0.5 

5590 

-984 

.0320 

0 

SAME  DATA  AS  SCENARIO 

ONE 

************  constant  G turn  (approximately) 

6.0 

1000 

-20 

1 .0000 

102 

6.5 

990 

-20 

1 . 005 

102 

1 

7.0 

980 

-19 

1.005 

102 

7.5 

970 

-19 

1.005 

102 

1 1 

( 

SCENARIO  THREE 

Time 

(seconds ) 

Range 
(meters ) 

Range  Rate 
(meters/Sec ) 

Angular  Rate 
( radi ans/Sec  ) 

Acceleration 

(P's) 

0.5 

1250 

-180 

.1920 

0 

1.0 

1165 

-154 

.2206 

0 

1.5 

1097 

-123 

.2495 

0 

1 j 

2.0 

1044 

-86 

.2752 

0 

i 

2.5 

1011 

-44 

.2934 

0 

3.0 

1000 

-1.8 

.3000 

9.2 

1 1 j 

3.5 

999 

-1.8 

.3004 

9.2 

4.0 

998 

-1.8 

.3007 

9.2 

4.5 

997 

-1.8 

.3009 

9.2 

1 j 
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experiences  the  same  deceleration  mentioned  in  scenario 
four.  This  scenario  is  illustrated  in  Figure  10.  The  time 
history  for  this  scenario  is  presented  in  Table  V. 

It  is  important  to  note  that  for  scenarios  two,  three, 
four  and  five  the  conversion  from  one  acceleration  profile 
to  another  was  modeled  as  an  instantaneous  change.  Although 
this  instantaneous  change  in  profiles  is  not  realistic,  it 
again  accentuates  the  response  characteristics  of  the 
system. 


Velocity  = 1000  m/sec 
Range  at  Crossover  = 1000  m 


0 t = 4.5  D = -\p CdMc 

CQ  = .25  W = 90  kg 

radius  = 15  cm  deceleration  = 7 G's 


Crossover 


t = 4.5 


Fig.  9.  Scenario  Four,  Flyby  to  Deceleration 


Velocity  = 1000  m/sec 
0 t = 2 sec  target  enters  15  G 


turn  toward  tracker 

0 rollout  target  decelerates  as 
in  scenario  four  , 


At  = 2 sec 


'Rollout 


Fig.  10.  Scenario  Five,  Flyby  to  Pursuit 
to  Deceleration 
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TABLE  IV 
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SCENARIO  FOUR 


Time 

(sec) 

Vel oci ty 
(m/sec ) 

Range 
(meters ) 

Range  Rate 
(m/sec ) 

Angular  Rate 
(rad/sec ) 

Accel erati on 

-(G,_s) 

0.5 

1000 

5590 

-984 

. 0320 

0 

SAME  DATA 

AS  SCENARIO 

ONE 

**********  decelerate  ********** 

4.5 

1000 

1803 

-832 

.3072 

-7.8 

5.0 

962 

1421 

-683 

.4760 

-7.2 

5.5 

927 

1135 

-439 

.7193 

-6.7 

6.0 

896 

1003 

-73.4 

.8896 

-6.3 

6.5 

866 

1062 

292 

.7675 

-5.9 

7.0 

838 

1271 

517 

.5190 

-5.5 

7.5 

812 

1559 

623 

.3339 

-5.2 

8.0 

787 

1884 

667 

.2219 

-4.9 

: 
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TABLE  V 


SCENARIO  FIVE 


T ime 
(sec ) 

Range 
(meters ) 

Range  Rate 
( in  / s e c ) 

Angular  Rate 
( rad/sec ) 

Velocity 
(m/sec  ) 

Acceleration 

(G's) 

0.5 

5590 

-984 

. 0320 

1000 

0 

1.0 

5099 

-980 

.0385 

1000 

0 

1.5 

4610 

-976 

.0471 

1000 

0 

*********  15.6 

G turn  toward  tracker  ********* 

2.0 

4123 

-970 

.0585 

1000 

15.6 

t 

2.5 

3635 

-981 

.0533 

1000 

15.6 

3.0 

3142 

-988 

.0454 

1000 

15.6 

3.5 

2646 

-996 

.0327 

1000 

15.6 

4.0 

2147 

-999 

.099 

1000 

15.6 

*********  ronout  and  decelerate  ********* 

4.5 

1652 

-973 

.0004 

973 

-7.4 

5.0 

1174 

-938 

.0009 

938 

-6.9 

5.5 

713 

-906 

.0008 

906 

-6.4 

i 
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V . Results  Analysis 


This  chapter  presents  the  basic  results  of  the  study 
with  an  analysis  of  what  the  results  signify.  Initially, 
an  analysis  of  the  parameter  Q33  versus  RMS  tracking  error 
is  presented  in  order  to  determine  the  optimum  value  of  Q ^ 3 * 
Next,  the  results  of  using  a computational  delay,  instan- 
taneous update,  and  prediction  update  in  generating  W ^ 
are  presented.  Subsequently,  with  Q 3 ^ set  at  the  optimum 
value  and  the  best  method  of  generating  W ^ used,  the 
results  indicating  that  Modification  Six  demonstrated  the 
best  performance  are  presented.  Finally,  a comparison  of 
the  results  from  the  faster  0.01  second  sample  rate  pseudo- 
range system  versus  the  slower  0.1  second  sample  rate  system 
is  presented. 

The  measurement  noise  in  this  simulation  was  modeled  as 
a white,  Gaussian,  zero  mean  sequence.  The  standard  devia- 
tion in  range  is  5 meters  and  in  pointing  angle,  10  micro- 
radians. The  same  sequence  was  used  for  each  simulation. 


Q33  Analysis 

Figure  11  is  a graph  of  Q33  versus  RMS  tracking  error. 
Curves  1 and  2 represent  the  time  averaged  RMS  tracking  per- 
formance durinq  1.5  seconds  of  the  100  G turn  in  scenario 
two,  with  and  without  noise  respectively.  From  these  curves 
a value  of  Q33  = 1.0  appears  to  be  optimal.  However,  curves 
3 and  4 represent  the  RMS  tracking  performance  during  the  non- 
accelerating  portion  of  scenario  two,  with  and  without  noise 


respectively.  These  curves  indicate  that  Q33  should  be  set 

as  low  as  possible  when  tracking  a non-accelerating  target. 

On  the  other  hand,  curves  5 and  6,  which  represent  the 

performance  during  1.5  seconds  of  the  10  G turn  in  scenario 

three,  indicate  a of  approximately  0,01  to  be  optimal. 

Since  actual  target  motion  is  unknown,  a tradeoff  must  be 

made  between  Q33  = 0.0  and  1.0,  because  to  choose  either 

one  exactly  would  causer  substantial  errors  if  the  target 

» 

motion  was  just  the  opposite.  From  analysis  of  Figure  11, 

Q33  = 0,1  was  chosen  as  the  optimum  tradeoff  value  to  handle 
both  highly  accelerating  and  non-accelerating  targets. 
Ideally,  the  value  of  Q33  should  be  set  according  to  the 
bandwidth  of  the  random  or  unknown  target  motion  as  indicated 
by  Fitts  (Ref  3:73). 


Generating  W 

The  results  of  the  discussion  in  Chapter  III  about 
the  various  techniques  of  generating  ^ are  shown  in 
Figures  12  through  26.  (Note:  Figures  12  through  67  are 
contained  in  Appendix  E and  are  plots  of  tracking  error 
versus  time.)  Table  VI  gives  a listing  of  the  RMS  and  peak 
tracking  errors  for  each  technique.  Originally,  Modifica- 
tion Two  was  used  for  this  analysis.  Later,  Modification 
Six  was  compared  to  deduce  what  conclusions  could  be  drawn 
for  it.  A close  observation  of  Figures  12  through  26  and 
of  Table  VI  yields  the  following: 
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TABLE  VI 

Performance  of  Techniques 
Used  to  Generate  M0 ^ ^ 


Modification  Peak  Error  RMS  Error  (microradians ) 


Two 

(microradians ) 1,5- 

6 sec  6-8  sec 

1.5  - 

* ★ * ★ 

Scenario 

One 

★ ★ ★ ★ 

Del ayed 

-35.5 

9.6 

13.4 

10.7 

Instantaneous 

20 

4.6 

9.1 

6.3 

Pred ict 1 on 

15.7 

4.4 

7.0 

5.3 

* * * * 

Scenario 

Two 

* if  it 

Delayed 

531  .9 

9.5 

296 

149.4 

Instantaneous 

533.9 

4.6 

297 

149.6 

Prediction 

527.0 

4.4 

294 

147.8 

★ ★ * fc 

Scenario 

Five 

Delayed 

96.2 

26.9 

Instantaneous 

82.0 

28.0 

Prediction 

95.4 

26.8 

Modi fi cat  ion 
S i x 


**** 

Scenario 

One 

★ *** 

Del ayed 

-32.8 

8.2 

12.8 

9.7 

Instantaneous 

15.1 

2.7 

6.5 

4.2 

Prediction 

15.3 

3.8 

6 . 6 

4.8 

**** 

Scenario 

Two 

Delayed 

388.9 

8.1 

207 

104.3 

Instantaneous 

390.4 

2.7 

208 

104.7 

Prediction 

406 

3.9 

215 

108.3 

40 
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1.  When  Modification  Two  was  used  aqainst  scen- 
ario one,  the  predicted  update  technique  showed 
better  overall  performance  in  terms  of  peak 

and  RMS  errors  than  the  computational  delay 
or  instantaneous  update  technique.  However, 
when  Modification  Six  was  used,  the  instan- 
taneous technique  demonstrated  the  best  per- 
formance . 

2.  When  Modification  Two  was  used  aqainst  scen- 
ario two,  the  predicted  update  technique  was 
again  superior.  On  the  other  hand.  Modifica- 
tion Six  showed  mixed  results.  The  computa- 
tional delay  technique  provided  slightly  lower 
peak  and  RMS  errors  durinq  the  10  G turn. 

Whereas,  the  instantaneous  technique  was 
substantially  better  in  the  unaccelerated 
portion  of  the  trajectory  than  the  other  two 
techniques. 

3.  Against  scenario  five,  the  results  are  incon- 
clusive. Figures  13  through  20  show  no  appreci- 
able difference.  For  this  reason.  Modification 
Six  was  not  compared  against  scenario  five. 

From  the  above  results,  it  appears  that  the  predicted 
update  technique  improves  the  tracking  performance  of  Modi- 
fication Two  over  that  obtained  using  the  "ideal"  instan- 
taneous update  technique.  However,  when  Modification  Six 
was  compared,  the  instantaneous  update  technique  provided 
the  best  overall  performance.  Throughout  the  remainder  of 
this  study  the  instantaneous  update  technique  was  used. 


Modification  Selection 

The  five  different  methods  of  generating  the  pseudo- 
range measurement  were  used  to  track  the  target  modeled  by 
scenario  two.  Figures  27  through  36  illustrate  the  tracking 
performance  attained  by  each  method  with  and  without  noisy 
measurements . 
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A comparison  of  Figures  28,  30,  32,  34,  and  36  show 
that  noise  has  relatively  the  same  effect  on  each  modifica- 
tion with  the  exception  of  Modification  Five  as  predicted 
in  Chapter  IV.  Noise  in  the  measurements  is  propagated 
more  by  the  differenced  range  rate  expression  than  by  any 
of  the  other  methods. 

A comparison  of  Figures  27  through  36  and  of  the 
tabulated  data  in  Table  VII,  clearly  shows  Modification 
Six  (using  the  estimated  range)  provides  the  best  overall 
tracking  performance. 

Performance  Com p a r i s o n 

This  section  presents  a compilation  of  the  most  signi- 
ficant results  of  this  study.  The  results  from  scenarios 
one,  two,  and  five  are  analyzed  here.  Table  VIII  provides 
a tabulation  of  peak  and  RMS  tracking  errors  used  in  making 
the  comparisons. 

Scenario  One.  Figure  37  illustrates  the  tracking 
error  obtained  when  a conventional  pointing  and  tracking 
system  is  used  to  track  target  motion  modeled  by  scenario 
one.  Peak  tracking  error  is  927  microradians  with  the  RMS 
tracking  error  being  492  microradians.  Figure  38  illustrates 
the  performance  obtained  when  the  0.1  second  Aided  Track 
system  is  used.  Peak  error  is  reduced  to  a maximum  value 
of  8.25  microradians.  RMS  tracking  error  is  also  reduced 
to  2.75  microradians. 
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TABLE  VII 


Performance 

of  Modifications 

Against  Scenario 

Two 

Modi f i cati on 

Peak  Error 
(microradians) 

RMS 

Flyby 

Error  (microradians) 
100  G Turn  Overall 

One 

533.0 

26.1 

286.6 

145.7 

One  with  Noise 

552.3 

36.9 

295.5 

151.6 

Two 

533.9 

4.6 

297.7 

149.6 

Two  with  Noise 

566.6 

11.3 

305.6 

153.8 

Three 

528.8 

4.3 

294.9 

148.3 

Three  with  Noise 

561.6 

11.1 

302.8 

152.4 

Five 

571.2 

4.4 

313.0 

157.5 

Five  with  Noise 

599.8 

17.8 

322.8 

162.8 

S i x 

390.4 

2.67 

208.4 

104.7 

Six  with  Noise 

385.5 

6.77 

207.8 

104.6 

Flyby  - 1.5  to  6 seconds 
100  G turn  - 6 to  8 seconds 
Overall  - 1.5  to  8 seconds 

Time  average  RMS  calculations  were  begun  1.0  seconds 
after  filter  was  activated  to  allow  for  transients  to  die 
out . 
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Results  shown  in  Figure  39  demonstrate  the  best 
possible  performance  attainable  by  the  system.  These 
were  obtained  using  Modification  Four.  The  peak  error 
is  increased  slightly  from  that  shown  in  Figure  38  to  a 
value  of  9.61  microradians.  The  RMS  tracking  error  is 
also  increased  to  3.3  microradians.  This  demonstrates  the 
degradation  in  performance  mentioned  by  Fitts  when  the 
sample  rate  is  increased  (Ref  1:91,97).  However,  the 
flyby  scenario  is  a constant  velocity,  zero  acceleration 
model.  The  filter  propagates  a constant  acceleration  over 
the  sampling  interval.  Therefore,  if  Q ^ 3 and  P33  are  set 
to  zero,  then  the  filter  model  matches  the  target  model. 
Figures  40  and  41  show  that  when  this  is  done,  the  degrada- 
tion is  no  longer  present.  On  the  other  hand,  the  faster 
system  shows  no  improvement  over  the  slower  system  other 
than  eliminating  the  estimation  oscillations. 

Figure  42  illustrates  the  performance  obtained  when 
Modification  Six  is  used.  Peak  tracking  error  is  increased 
even  more  to  15.2  microradians  and  the  RMS  tracking  error 
to  4.23  microradians.  This  increased  degradation  is  duo 
to  the  use  of  the  pseudo- range  measurements. 

Scenario  Two.  Figures  43  through  48  illustrate 
the  results  obtained  when  target  motion  was  modeled  by 
scenario  two.  A close  look  at  Figure  43  clearly  shows  that 
a conventional  unaided  pointing  and  tracking  system  produces 
essentially  zero  tracking  error  during  the  100  G,  constant 
radius,  constant  velocity  turn.  On  the  other  hand.  Figures 


44  through  46  demonstrate  the  fact  that  any  aiding  process 
results  in  relatively  large  errors  during  this  same  100  G 
turn.  This  is  due  to  the  basic  filter  model  assuming  con- 
stant acceleration  over  the  sample  interval,  which  is  not 
true  in  this  case.  However,  aiding  practically  eliminates 
the  large  tracking  error  prior  to  the  100  G turn  as  seen  in 
Figure  43. 

A comparison  of  Figure  44  and  45  shows  that  a definite 
improvement  in  the  tracking  performance  of  the  0.1  second 
system  during  the  100  G turn  is  possible  by  increasing  the 
sampling  rate  of  the  system.  Figure  45  however,  is  a graph 
of  the  tracking  performance  of  the  optimal  modification 
using  actual  laser  rangefinder  measurements  every  0.01 
seconds.  A comparison  of  Figure  44  and  46  illustrates  that 
the  use  of  pseudo-range  measurements  to  increase  the  sampl- 
ing rate  still  results  in  an  improvement.  The  peak  tracking 
error  during  the  100  G turn  is  reduced  from  699  microradians 
in  Figure  44  to  390  microradians  in  Figure  46.  Similarly, 
the  RMS  tracking  error  during  the  100  G turn  is  reduced 
from  369  microradians  for  the  unmodified  case  to  208 
microradians  for  the  system  using  Modification  Six. 

Figures  47  and  48  show  the  effect  noisy  measurements 
have  on  the  0.1  second  system  and  the  Modification  Six 
system  respectively.  Roth  appear  to  be  affected  to  about 
the  same  degree. 

Another  interesting  result  can  be  seen  by  comparing 
the  RMS  tracking  errors  for  the  non-accelerating  portion 
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(that  portion  prior  to  crossover)  of  scenario  two  shown  in 
Table  VIII.  The  0.1  second  Aided  Track  system  has  an  RMS 
error  of  2.04  microradians  for  this  portion  of  the  trajec- 
tory; whereas,  the  0.01  second  Modification  Four  system  has 
only  a 1.88  microradian  RMS  error.  This  is  in  contradiction 
to  the  degradation  prediction  found  when  the  entire  trajec- 
tory of  scenario  one  was  considered.  However,  Modification 
Four  uses  actual  laser  rangefinder  measurements.  When  the 
pseudo-range  measurements  in  Modification  Six  are  used,  the 
RMS  tracking  error  increases  up  to  2.67  microradians.  The 
primary  differences  between  Modification  Four,  Modification 
Six,  and  the  unmodified  systems  can  be  seen  by  comparing 
Figures  38,  39,  and  41  from  scenario  one  and  Figures  44, 

45,  and  46  from  scenario  two . Modification  Four  totally 
eliminates  the  range  correction  oscillations  in  the  track- 
ing error  curves  seen  in  the  unmodified  system  and  the 
Modification  Six  system. 

Scenario  Five.  The  results  obtained  when  target 
motion  was  modeled  by  the  more  realistic  scenario  five  are 
depicted  in  Figures  49  through  51.  The  tracking  performance 
of  a conventional  unaided  system  is  shown  in  Figure  49. 

The  jumps  in  the  curve  at  2.0  seconds  and  4.4  seconds  are 
from  the  16  G turn  toward  the  tracker  and  the  deceleration 
model  respectively.  The  small  jump  at  5.0  seconds  is  due 
to  the  target  reentering  the  16  G turn  profile  momentarily 
because  of  a simulation  modeling  error.  The  peak  tracking 
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Scenario  One 

Peak  error 
(microradians ) 

RMS  error 
(microradians ) 

No  aiding 

927 

492 

0.1  second  rate 

8.25 

2.75 

Modification  Four 

9.61 

3.3 

Modification  Two 

20 

6.34 

Modification  Six 

15.2 

4.23 

Scenario  Two 

Peak  error 

RMS  error 

in  flyby  in  turn 

No  aiding 

927  (in  flyby) 

365 

20.9 

0.1  second  rate 

699  (in  turn) 

2.04 

369 

with  noise 

770  (in  turn) 

8.08 

366 

Modification  Four 

356  (in  turn) 

1.88 

277 

Modi f i ca t i on  Two 

533.9  (in  turn) 

4.63 

297 

Modification  Six 

390.4 

2.67 

208.4 

with  noise 

385.5 

6.77 

207.8 

Scenario  Five 

Peak  error 

RMS  error 

No  aiding 

-100 

37 

0.1  second  rate 

240 

67 

Modification  Two 

82 

•28 

Modification  Six 

95 

27 

error  in  this  case  is  approximately  -100  microradians , 
which  occurs  just  as  the  target  is  rolling  out  aimed  at 
the  tracker.  Once  the  target  rolls  out  and  begins  to 
decelerate,  the  conventional  system  appears  to  track  very 
accurately.  The  overall  RMS  tracking  error  is  37  micro- 
radians  . 

The  0.1  second  Aided  Track  system  results  are  shown 
in  Figure  50.  This  system  appears  to  have  great  diffi- 
culty at  close  ranges  and  head-on  attacks.  The  peak 
tracking  error  occurs  during  the  head-on  attack  portion  of 
the  trajectory,  the  value  being  approximately  240  micro- 
radians.  The  overall  RMS  tracking  error  is  67  microradians, 
which  is  worse  than  the  conventional  system. 

The  performance  obtained  by  using  Modification  Six  is 
shown  in  Figure  51.  A substantial  improvement  over  the 
conventional  system  and  the  0.1  second  Aided  Track  system 
can  be  seen.  The  peak  tracking  error  is  only  95  micro- 
radians and  the  RMS  error  is  only  27  microradians.  This 
system  tracks  the  head-on  attack  much  better  than  the  0.1 
second  aided  system  and  also  handles  the  16  6 turn  better 
than  either  the  conventional  system  or  the  0.1  second  Aided 
Track  system. 

The  remaining  Figures  52  through  67  are  provided  as 
supportive  material . 
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Cone  1 us i ons 


1.  A value  of  Q33  = 0.1  provided  the  optimum  trade- 
off at  the  0.01  second  sample  rate.  This  conclusion  is 
consistent  with  that  made  by  Fitts. 

2.  Generating  Wfl ^ ^ by  the  instantaneous  update 
technique  is  the  "ideal"  situation  and  provided  the  best 
performance  of  the  three  computational  techniques  when 
Modification  Six  was  used.  However,  the  predicted  update 
technique  performed  almost  as  well.  The  instantaneous 
technique  is  not  realisitic  in  terms  of  hardware  whereas 
the  prediction  technique  is. 

3.  Modification  Six,  which  uses  the  estimate  of  range 
from  the  filter  as  a pseudo-range  measurement,  demonstrated 
the  best  performance  of  all  the  modifications  tested. 

4.  The  faster  sample  rate  results  in  a degradation  in 
performance  when  tracking  simple  target  motion,  such  as  the 
flyby  scenario,  with  Q ^ 3 = 0.1  and  P33  = 1.0.  However, 
when  Q33  = P33  = 0.0,  there  is  no  degradation  but  there  is 
no  improvement  either. 

5.  The  faster  sample  rate  does  result  in  an  improve- 
ment in  tracking  performance  when  tracking  highly  maneuver- 
able  targets,  especially  when  the  target  is  attacking  the 
tracker. 


Recommon da ti ons 

1.  From  the  Q33  analysis,  a tuneable  Q33  is  desirable, 
dependent  upon  the  maneuverability  of  the  target  being  tracked. 
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2.  The  predicted  update  technique  for  generating  W 
should  be  incorporated  in  the  system. 

3.  When  tracking  flyby  or  non-accelerating  type  tar- 
gets, the  faster  sample  rate  system  should  not  be  used. 

4.  When  tracking  highly  maneuverable  targets,  Modifi- 
cation Six  should  be  incorporated  to  improve  the  performance. 

5.  In  order  to  implement  this  algorithm,  the  follow- 
ing steps  should  be  taken. 


a .  ) Set  AT  = 0. 01 

b .  ) Set 


H = 


0 0 0 

0 1 0 


between  laser  rangefinder  measurements  and 


H = 


1 


0 O' 


[0  1 OJ 

at  laser  range  measurement  time. 

Remark:  The  Kalman  gains  can  still  be  precomputed  and 

stored,  since  the  ^matrix  is  changed  in  a 
periodic  manner. 

c. )  Use  the  estimate  of  range  from  the  filter  as 

the  pseudo-range  measurement  between  the 
laser  range  updates. 

d. )  Use  the  predicted  update  technique  for 

generating  W.^  in  order  to  improve  the 
performance  of  the  system. 


50 


Bibliography 


i 

* 


1.  Fitts,  J.  M.  HAT  Aided  Tracking  Analysis.  Report  No. 

P71-2154.  "Culver  City,  California:  Hughes  Aircraft 
Company,  June  1971. 


i 


f ' 


2.  Fitts,  J.  M.  A Strapdown  Approach  to  HAT  Aided  Tracking. 

Report  No.  >71-351.  Culver  City,  California: 

Hughes  Aircraft  Company,  September  1971. 

3 . Fitts,  J . M . H AT  Position  and  Rate  Aided  Tracking 

Analysis.  Report  No.  P71-471.  Culver  City, 
California:  Hughes  Aircraft  Company,  December 

1971  . 


4. 


5. 


Fitts,  J.  M.  "Aided  Tracking  as  Applied  to  High  Accuracy 
Pointing  Sy s t em s . " IEEE  Transactions  on  Aerospace 
and  Electronic  Sy stems  . 9:35  0-368  (May  1 973  ) . 


User  Manual  for  QSUAPT-3.0  and  FROAPT,  Report  ER-73-R- 
22-003.  Stillwater,  Oklahoma:  Oklahoma  State 
University,  January  1974. 


Supplemental  Bibliography 


Air b or n e_  Po i n t i no  and  Tracking  System  (U),  AFWL-TR-72- 
~124,  Kirtland  Air  Force  Base,  New  Mexico:  Air 
Force  Weapons  Laboratory:  1972. 

Negro,  0.  E.  and  Robert  R.  Wilde,  "An  Aid  for  Improving 
The  Precision  Pointing  of  the  FTT."  Laser  Digest, 
AFWL-TR-72-243 , Kirtland  Air  Force  Base,  New  Mexico: 
Air  Force  Weapons  Laboratory,  December  1972. 

Sebesta,  H.  R.  and  L.  R.  Ebbesen.  Studies  of  Advanced 
Simulation  Concepts.  Stillwater,  Oklahoma: 

Oklahoma  State  University,  October  1974. 


■I. 


Appendix  A 


Program  Listing:  A Rate  Aided  Pointing  and  Trackim 
Algorithm  for  a Two-Dimensional  System. 
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Truth  Model  Formulation 


In  order  to  model  the  target  motion  in  two  dimensions, 
all  that  is  needed  is  the  target  acceleration  in  the  Y^, 

Y^  coordinate  system  defined  in  Chapter  II,  and  a set  of 
initial  conditions.  The  differential  equations  can  then 
be  written  in  matrix  form  as 
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where  A^  and  A ^ are  the  modeled  target  acceleration  expres- 
sions in  the  Y^  , Y^  directions  respectively. 

The  differential  equations  describing  the  gimbal 
dynamics  can  be  written  by  inspection  of  the  closed  loop 
servo  portion  of  Figure  3.  By  looking  at  the  C^s)  and 


G^($J  transfer  functions  the  following  matrix  equation  is 


obtained: 


1 

0 


-w 
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1 

■?r  W. 


n„ 
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„P 
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w 't 

L E 

n 

cmd 


IB3) 


6A 
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r 


r 


wmm  t =t:t 


~~ - ■ - « 


And  from  the  Gj(S)  transfer  function  the  following  propor- 
tional plus  integral  compensation  relationship  is  obtained 


W,  = (700/18  + 7 0 0 / S ) r 


(84) 


Now  let 


then 


but 


C = ( 700/S  )t 


C = 70 Or  = 700  (nT  - np) 


n j = Arctan  ( Y ^ / V i ) 


(B5) 


(B6J 


(B7) 


Therefore , 


WA  = 700  (Arctan  (Y2/Y])  - n )/ 1 8 + C 


(B8) 


By  combining  Eq  (Bl),  (B2),(B6)»  and  (B8)  with  the  relation- 
shi  p 


Wcmd  = WA  + Waid 


(89) 


the  overall  truth  model  matrix  can  be  written  as  seen  in 

Eq  (8). 
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Appendix  C 


i 


\ 


Steady-State  Perturbation  Ana  1 ys i s 

Ihe  use  of  pseudo-range  measurements  to  increase  the 
sample  rate  of  the  Aided  Track  algorithm  is  susceptible  to 
the  problem  of  error  propagation.  Modification  Two  is  the 
worst  case  since  estimates  of  range  rate  are  used  every 
0.01  seconds  to  propagate  the  last  range  measurement.  This 
case  is  analyzed  here . 

Range  is  propagated  by  the  equation 

R' (K+l ) = R' (K)  + ft(K)AT  (Cl ) 

where  the  prime  denotes  pseudo-range  and  R'(0)  is  the  last 
laser  range  measurement.  Assuming  that  the  laser  range 
measurement  is  perfect,  the  estimate  of  range  rate  R and 
the  pseudo-range  terms  can  be  expressed  and  the  sum  of  their 
true  values  plus  some  error. 


I I I 


R(K)  = 

R(K)  = 6R(K),  K = 0 , 1 ....  ,9 

(C?) 

R'U)  = 

R(K)  + <$R  1 (K)  , K=1  ,2, . . . ,9 

( C 3 ) 

where  6R(0)  = 0. 

Eq  (Cl)  can  be  rewritten  as 

I I 

R ' ( K + l ) = R ( K ) + 6R(K)  + R{ K) AT  + sR(K)AT,  K=l,..,9  (C4) 

But  the  first  and  third  terms  represent  the  true  range. 
Therefore,  the  error  in  the  pseudo-range  propagation  can 
be  written  as 


70 


* i 


6R* 


(K+l ) = «R' (K)  4 a^(K)AT,  K=1 9 


IC5) 


were  gR ' (0)  = 0. 

The  pseudo-range  measurements  are  used  in  Eq  (6) 
through  (9)  which  are  restated  here. 


X,  ■ R'  cos  n 
111 


Xp  = R'  sin  n 
^ in 


(C6) 


X,  = R cos  n - R'  W sin  nm 
in 


Xp  = "R  sin  n + R‘  W cos  nm 
m 


Assuming  that  n and  W are  known  perfectly,  the  perturba- 
tions caused  in  Eq  ( C 6 ) by  the  use  of  pseudo-range  measure- 
ments are  expressed  as 


«X,  = C «R  • 

m 

+ 6$AT] 

cos 

nm 

6Xp  = E6R 1 
m 

1 + fi^AT] 

s i n 

nm 

6X,  = R 

m 

cos  nm  - 

C 6 R • 

+ 5^ at!  W sin  n 

«X2  = R 

sin  nm  + 

[ 6 R 1 

+ s'RatI  W cos  n 

IC7) 


Assume  the  nominal,  steady-state  conditions  to  be 


nm  * 90° 


W = 1 radian/second 
AT  = 0.01  seconds 
fi  t\  (0)  = 0.01  meters/second 
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Then  Eq  (C7)  becomes 


6 X,  = 0 
m 


6 X 2 = 6 R ' + 


6 X -I  =-[6R‘  + S'RatJ 
m 


dX2  = 6^ 
m 


( C8 ) 


Eq  ( C8 ) is  further  simplified  if  only  the  errors  after  the 
first  iteration  are  considered  (i.e.  6 R 1 (0)  = 0). 


6 X -j  = 0 


6X2  = 0.001 


6 X -j  = -0.0001 


6X2  = 0.01 


( C9 ) 


These  perturbations  are  processed  through  the  filter  by 
the  following  set  of  equations: 


6X,+  = <5X,~  + G,  , (6  X,  - 6 X 
1 1 m 


1")  + G^U^  -6^") 

6X?+  = 6X2‘  + G21(6X2  - «X2")  + G22  U X2  -6^2") 

L c £m  m 

6^1  + = 6$,  + G-,,(6X-j  - 6X-,  ) + G 2 ( 6 X -j  -6^-|  ) 

m m 

6^o  = G?  •,  (6  X2  - 6X2  ) + G22  (6  X2  -6  ^2  ) 

m m 


(CIO) 
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Substituting  in  the  steady-state  Kalman  gains,  the  values 

expressed  in  Eq  (C9),  and  assuming  6^(0)  = 6X2(0)  = 

« • 

6X^(0)  = 6X2(0)  = 0 Eq  (CIO)  becomes 

6X]  = ( . 26 ) ( 0 ) + ( .03) (-.0001  ) = 3.0E-6 

6X2  = ( . 26 ) ( . 0001 ) + ( . 03) ( . 01 ) = 3.26E-4 

A (C11) 

6X]  = (.48)(0)  + (.11)1-. 0001)  = 1.1E-5 

6X2  = ( . 48)  ( . 0001  ) + ( . 1 1 ) ( . 0 1 ) = 1.148E-3 

The  perturbati ons  expressed  in  Eq  ( C 1 1 ) are  trans- 
formed into  the  change  of  the  error  in  the  estimate  of 
range  rate  by 

A 

A6R  = 1 1 

(6X1  + 6X2  7 (C12) 

= 1 .148F-3 

Therefore,  the  new  error  in. the  estimate  of  range  rate 
becomes 

I 

(1  ) = <5^(0)  + A6R(0) 

= .01  + .001  148  ( Cl  3 ) 

= .01148 

Thus  it  has  been  shown  that  the  use  of  pseudo-range 


measurements  creates  a stability  problem  of  errors  in  the 
estimates  of  range  rate  generating  larger  errors  in  the  next 
estimate  of  range  rate. 
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Program  Listings  for  Each  Modification  Method  Tested 
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